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Motivated by recent studies on the quasi-one-dimensional (ID) antiferromagnet BaCo2V208, 
we investigate the Tomonaga-Luttinger-Iiquid properties of the 10 5 = 1/2 XX Z Heisenberg- 
Ising model in magnetic fields. By using the Bethe ansatz solution, thermodynamic quantities 
and the divergence exponent of the NMR relaxation rate are calculated. We observe a magne- 
tization minimum as a function of temperature (T) close to the critical field. As the magnetic 
field approaches the critical field, the minimum temperature asymptotically approaches the 
universal relation in agreement with the recent results by Maeda et al. We observe T-linear 
specific heat below the temperature of the magnetization minimum. The field dependence of 
its coefficient agrees with the results based on conformal field theory. The field dependence of 
the divergence exponent of the NMR relaxation rate with decreasing temperature is obtained, 
indicating a change in the critical properties. The results are discussed in connection with 
experiments on BaCo2V2 08. 
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1. Introduction 

One-dimensional (ID) quantum spin systems have 
been investigated for many years. Their fascinating as- 
pects have been revealed by theoretical and experimen- 
tal studies. In particular, ID gapped spin systems have 
attracted a large amount of attention.^ There are vari- 
ous gapped spin systems including Haldane-gap systems, 
two-leg ladders, and bond-alternating systems. It was 
argued that the characteristics of each system appear 
conspicuously in their critical and dynamical properties 
above the critical magnetic field (He) where the gap is 
closed and low-energy properties can be described as a 
Tomonaga-Luttinger liquid (TLL).^ Such characteristic 
behavior depending on the model can be observed in, 
for example, the NMR relaxation rate.^'^ In fact, the 
divergent behavior of the NMR relaxation rate with de- 
creasing temperature was observed for S — 1 Haldane- 
gap material (CH3)4NNi(N02)3'*'^ and 5=1/2 bond- 
alternating chain material pentafluorophcnyl nitronyl ni- 
troxide (FsPNN).^ The field dependence of the diver- 
gence exponent was discussed in comparison with theo- 
retical results.''"^ 

Recent experimental advances have allowed the study 
of the field dependence of thermodynamic quanti- 
ties in ID gapped spin systems in precise com- 
parison with the theory. The specific heat of the 
quasi-lD S — 1 bond-alternating antiferromagnet 
Ni(C9H24N4)(N02)C104 was measured in magnetic 
fields. ^° It was found that the low-temperature specific 
heat is proportional to temperature (T), indicative of 
the gapless dispersion relation. The field dependence of 
its coefficient was compared with the result based on 
conformal field theory (CFT).^^'^^ Good agreement was 
obtained. Other noticeable features of the TLL observed 
in its thermodynamic quantities were pointed out quite 
recently. The minimum of the magnetization appears as 



a function of temperature close to H = He in ID gapped 
spin systems. ^'^ The magnetization minimmn marks the 
important temperature below which the TLL is pre- 
served. It was demonstrated further that the tempera- 
ture of the minimum approaches the universal relation 
as H approaches He from above. It has been difficult to 
determine the temperature region for the TLL from the 
model Hamiltonian itself. The magnetization minimum 
plays an important role in determining the temperature 
region of the TLL not only in calculations but also in 
experiments. 

The quasi- ID antiferromagnet BaCo2V208 has been 
investigated intensively in recent years. From 
the magnetization measurement, it was shown that 
BaCo2V208 is well described by the ID S' = 1/2 XXZ 
Heisenberg-Ising model. The result has stimulated fur- 
ther investigations on the novel properties caused by the 
TLL. Since the Neel-like ground state of the model can be 
regarded as a charge-density-wave (CDW) insulator with 
quasi-long-range order from the viewpoint of the spinless 
fermion, it is expected that the dominant spin correlation 
changes from longitudinal incommensurate to transverse 
staggered with increasing magnetic field when He < H . 
The change in the critical property can be detected via 
the NMR relaxation rate.^'^'^ It was argued previously 
that the change in the dominant spin correlation was de- 
tected for F5PNN via the NMR relaxation rate.^ How- 
ever, the conclusion is not yet settled. BaCo2V208 is a 
potential candidate for observing the change in the dom- 
inant spin correlation in magnetic fields. 

In this paper, we investigate the TLL properties of 
the ID S = 1/2 XXZ Heisenberg-Ising model in mag- 
netic fields by using the Bethe ansatz solution. Thermo- 
dynamic quantities such as the magnetization, suscepti- 
bility, and specific heat are calculated. The temperature 
dependence of the magnetization is first calculated. From 
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the temperature of the magnetization minimum, we eval- 
uate the temperature region of the TLL in a given mag- 
netic field close to H = He- The results are used for dis- 
cussing the TLL features obtained in other calculations. 
The field dependence of the coefficient of the T-linear 
specific heat is investigated. The results are compared 
with the field dependence of the coefficient obtained from 
the CFT relation 7r/3u with v being the excitation veloc- 
ity. Furthermore, the divergence exponent of the NMR 
relaxation rate is calculated as a function of the magnetic 
field. In §2, we outline the model and the method. In §3, 
we present the numerical results, which are discussed in 
connection with the thermodynamic and critical proper- 
ties of BaCo2V2 08. Section 4 is a summary. 

2. Model and Method 

We consider the ID S = 1/2 XXZ Hcisenberg-Ising 
model in a magnetic field 

i 

- giXBHY^St. (1) 

i 

where J > 0. We set g^-Q — 2. By the Bethe ansatz 
solution, the critical field {He) and the saturation field 
[Hs) are given as^® He = JK{u')u' smh(f)/(f) and Hg = 
J(l -I- A)/2, respectively, where K{u') is the complete 
elliptic integral with its modulus u', and A — cosh0 
(0 > 0). The modulus u' is determined by K{u')/ K{u) — 
with u' — ^/l — V?. The measured magnetization 
of BaCo2V208 was fitted to the magnetization curve 
at T = obtained from the Bethe ansatz solution, 
yielding A = 2.17 and J = 30.0 K (/cb = 1).^° In this 
case. He ~ 0.25J and Hg ^ 1.58J. It was pointed out 
that the ferromagnetic next-nearest-neighbor interaction 
J' = — 0.065J is required to achieve quantitatively pre- 
cise agreement with the results of ESR measurements.^^ 
However, this correctional interaction does not act as 
frustration in the Neel-like ground state of the Hamil- 
tonian in eq. (1). In the present calculations, therefore, 
we neglect its small effect. In the following, J is used in 
units of energy. 

The formalism for the Bethe ansatz solution at a fi- 
nite temperature was derived by Gaudin.^^ The basic 
equations for equilibrium thermodynamics are written 
as an infinite set of nonlinear integral equations for the 
pseudoenergy. To obtain accurate numerical results, we 
employ the technique developed for the impurity Ander- 
son model, where thermodynamic quantities were cal- 
culated without numerical differentiation of the free en- 
ergy functional. To accomplish the numerical calculation, 
we necessarily cut off an infinite set of the nonlinear in- 
tegral equations at a finite number, which is set to be 35. 
The iterative calculations for the pseudoenergies are car- 
ried out until they all converge within a relative error of 
10"^. Using the obtained pseudoenergies, their deriva- 
tives are calculated in the same way within a relative 
error of lO"**. 

Using the thus-obtained pseudoenergies and their 
derivatives with respect to temperature, we calculate the 



temperature dependence of thermodynamic quantities in 
magnetic fields. 

3. Results and Discussion 

3.1 Magnetization minimum and susceptibility 
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Fig. 1. (Color online) The magnetization in a given magnetic field 
as a function of temperature. 
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Fig. 2. (Color online) The temperature of the magnetization min- 
imum as a function of the magnetic field. The solid circles rep- 
resent Tm and the solid line expresses the universal relation, 
Tni = xo{H — He), [xQ ^ 0.76238). Inset: the velocity as a func- 
tion of the magnetic field. He ~ 0.25 and Hb ~ 1.58. 



We show the magnetization as a function of tempera- 
ture in Fig. 1; the magnetization has a minimum close to 
H = He- The decrease in m with increasing temperature 
is caused by the positive dv/dH in the low-temperature 
expansion for the magnetization based on CFT.^"^ By us- 
ing the Bethe ansatz solution at T=0, the velocity is ob- 
tained as V = [de{X)/d\]/[2np{X)]\x=B, where £(A) and 
p(A) are the dressed energy and the distribution func- 
tion for the rapidity A obtained by the following integral 
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equations 

e{X) - £o(A) 

p(A)-ai(A) 



dX'a2i\- \')e{X'), (2) 



dA'a2(A- A')p(A'), (3) 



-B 



with eo(A) = 2H — (27rsinh(/)/(/))ai(A) and a„(A) = 
sinh(n0)/[cosh(n(/)) — cos A]. The cutoff B is obtained by 
e{B) — 0. As shown in the inset of Fig. 2, the velocity in- 
creases immediately just above H — He and takes a max- 
imum around H = 0.63. The magnetization minimum 
appears when He < H < 0.63 in agreement with the the- 
oretical analysis. The temperature of the magnetization 
minimum (Tm) vs the magnetic field is plotted by the 
solid circles in Fig. 2. As the magnetic field approaches 
the critical field, Tm asymptotically approaches the uni- 
versal relation, = xo{H-He), {xo ~ 0.76238). Since 
the TLL comes into existence below Tm, the result helps 
us to discuss the temperature region for the characteris- 
tic TLL behavior seen in thermodynamic quantities and 
critical exponents. 




respectively. Both numerical and experimental results es- 
sentially agree with each other. 

3.2 Specific heat 
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Fig. 4. (Color online) The temperature dependence of the specific 
heat in magnetic fields. The lower panel represents the specific 
heat in a narrow range. 



Fig. 3. (Color online) The temperature dependence of the sus- 
ceptibilities below {H = 0.05), close to {H = 0.3), and above 
(H = 0.6) the critical field (H^ ~ 0.25). 



The susceptibilities in magnetic fields below, close to, 
and above the critical field He ^ 0.25 are shown in Fig. 3. 
The susceptibility for H — 0.6 approaches a finite value 
at T = with decreasing temperature, while the suscep- 
tibility for H — 0.05 approaches zero, reflecting the gap. 
We compare the numerical results with the experimen- 
tal ones. The susceptibility was measured in magnetic 
fields of H = LOT, H ^ 4.4 T, and H = 9.0 T.^^ Com- 
paring the critical field {He ~ 0.25) and the saturation 
field {Hs ~ 1.58) in our calculation with the observed 
values, He ^ 3.9 T and H^ ^ 22.7 T,!^ we evaluate that 
the three magnetic fields in the experiment correspond 
to H ^ 0.05, 0.30, and 0.60, respectively, in our model. 
In the experiments, the maxima of the susceptibilities 
were observed around T = 30 — 40 K in the three mag- 
netic fields. In Fig. 3, the maxima appear at T ~ 1.1 
for H = 0.05, T - 1.2 for H = 0.3, and T - 1.3 for 
H = 0.6, which correspond to T = 33 K, 36 K, and 42 K, 
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Fig. 5. (Color online) The coefhcient of the T-linear specific heat 
in magnetic fields. The solid circles represent the results obtained 
numerically and the solid line represents the result based on the 
CFT relation. 



The specific heat for H > He is shown in Fig. 4. 
When a quantum phase transition from a gapped state 
to a gapless state occurs in a ID quantum many-body 
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system, the shoulder or small peak structure may ap- 
pear for the low-temperature specific heat close to the 
quantum critical point. This behavior was first observed 
for the ID Hubbard model close to half-filling.^'^' Also 
for S = 1/2 two-leg spin ladders, the low-temperature 
peak of the specific heat was obtained for H > H^.^^ 
Ill the present system, wc; observe a shoulder structure 
at low temperatures for H = 0.35 and 0.6. It is consid- 
ered that the shoulder structure is caused by the con- 
tributions of both the T-linear specific heat at low tem- 
peratures, which is characteristic of the TLL, and the 
Schottky-type specific heat. Indeed, the T-linear specific 
heat emerges below T ~ and the shoulder structure 
appears slightly above it: Tm = 0.05 for H = 0.35 and 
Tm = 0.12 for H = 0.6. In Fig. 5, we compare the coef- 
ficient of the T-linear specific heat obtained by the low- 
temperature calculations with that obtained by the CFT 
relation 7 = 7r/3u.^^'^^ The former is represented by the 
solid circles and the latter is represented by the solid line. 
Both field dependences agree with each other. We have 
confirmed quantitatively that the shoulder structure of 
the low-temperature specific heat is caused by the TLL 
nature. 

Note that the low-temperature specific heat at H = He 
is proportional to VT,^^ reflecting the wave- number 
(q) dependence of the lower boundary of the excita- 
tion continuum: e{q) = —H^ -\- H^^ 1 + [u/u'Y sin^ q ~ 
{Hc/2){u/u'Yq^ for g ~ 0. Unfortunately, we could not 
reach such a low temperature region as to be able to 
find \/T dependence. It is a challenging issue to obtain 
C ~ \/T in a low temperature region at H = He- 

We find that the shoulder structures appear at T < 
0.15 for H = 0.35 and at T < 0.2 for H = 0.6. The mag- 
netic fields H = 0.35 and 0.6 are evaluated as H = 5.3 T 
and H = 8.8 T, respectively. Using J = 30.0 K, the 
temperature regions for the shoulder structure are eval- 
uated as T < 4.5K for H = 5.3T and T < 6K for 
H = 8.8 T. According to experiments on the specific 
heat for BaCo2V208,^^ no phase transition to the or- 
dered state takes place at least at T > 1 K for iJ > 5 T. 
Therefore, the T-linear specific heat as well as the shoul- 
der structure are expected to be observed. By comparing 
the observed field dependence of 7 and the results shown 
in Fig. 5, the TLL behavior of BaCo2V208 can be dis- 
cussed quantitatively. 

3.3 NMR relaxation rate 

We next investigate critical properties for < H < 
Hg. The critical exponent of the longitudinal spin cor- 
relation function is given by r]^ = 2[Z{B)]'^, where 

Z{X) = 1 - (l/27r) d\'a2{\ - \')Z{\').'^^ The criti- 
cal exponent 77^ takes the value of 1/2 at = He- As 
the magnetic field increases, rj^ increases monotonically 
and takes the value 2 at H = Hs- From the viewpoint of 
the critical properties, the gapped state for H < can 
be regarded as the CDW insulator with quasi-long-range 
order, while the fully spin-polarized state for Hg < H 
can be regarded as the band insulator. Since the weak 
universality 7/"^ • ry^ = 1 is preserved in the TLL with t]^ 
being the critical exponent of the transverse spin corre- 
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Fig. 6. The field dependence of the divergence exponent of the 
NMR relaxation rate. The V-like behavior represents the change 
in the critical property in magnetic fields. 



lation function, the dominant spin correlation changes at 
r?^ = 1 for iJc < < i?s-^'^'^ In the TLL of ID gapped 
spin systems in magnetic fields, the NMR relaxation rate 
(1/Ti) shows divergent behavior, which takes the form of 
1/Ti = A{H)T~°'. The divergence exponent a depends 
on the magnetic field and is described in terms of the 
critical exponent of the dominant spin correlation func- 
tion, reflecting the characteristic critical properties. The 
coefficient A{H) depends not on temperature but only 
on magnetic field. 

We find that 77^ = 1 at iJ = 0.82, which can be c;val- 
uated as H = 12.0 T. Accordingly, the incommensurate 
spin-density correlation is dominant with a = 1 — r]^ 
for 3.9T < H < 12.0 T, while the transverse stag- 
gered spin correlation is dominant with a = 1 — r]^ for 
12.0 T < H < 22.7 T. At if = 12.0 T, the divergent 
behavior of the NMR relaxation rate vanishes, because 
= jy^ = 1. In Fig. 6, we show the field dependence of a. 
The temperature region for the TLL around H = 0.82 
(12.0 T) can be roughly estimated as T < 0.13 (3.9 K) 
from the extrapolation of the results shown in Fig. 2. In 
the experiments on the specific heat, the phase transition 
to the ordered state was not observed until T = 0.4 K 
for H > lOT.i'^ Therefore, it is expected that the V- 
like behavior of a caused by the change in the critical 
property is observed below 3.9 K around H = 12.0 T for 
BaCoaVaOg. 

4. Summary 

We have investigated the TLL properties of the ID 
S = 1/2 XX Z Heisenberg-Ising model in magnetic 
fields. The thermodynamic quantities and the divergence 
exponent of the NMR relaxation rate have been calcu- 
lated by using the Bethe ansatz solution. We have found 
a magnetization minimum close to H ~ He. As H ap- 
proaches H = He the minimum temperature Tm, below 
which the TLL comes into existence, approaches the uni- 
versal relation asymptotically. We have observed the T- 
linear specific heat for T < T^. The field dependence of 
its coefficient agrees with the results based on the CFT 
relation. We have also demonstrated the V-like field de- 
pendence of the divergence exponent of 1/Ti, which re- 
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fleets the change in the critical property. The temper- 
ature and magnetic field regions to observe the above 
features in BaCo2V208 have been evaluated. 
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